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£S) ■ Abstract. The counts-in-cells and the two-point angular correlation function method 

t-H ' are used to test the randomnesses in the angular distributions of both the all gamma- 

ON . ray bursts collected at BATSE Catalog, and also their three subclasses ("short", "in- 

termediate", "long"). The methods eliminate the non-zero sky-exposure function of 
• BATSE instrument. Both tests suggest intrinsic non-randomnesses for the intermedi- 

ate subclass; for the remaining three cases only the correlation function method. The 
confidence levels are between 95% and 99.9%. Separating the GRBs into two parts 
("dim half and "bright half, respectively) we obtain the result that the "dim" half 
CZ2 ■ shows a non-randomness on the 99.3% confidence level from the counts-in-cells test. 
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>< : INTRODUCTION 

.«?: 

Recently, two different articles [1,2] simultaneously suggest that the earlier se- 
paration of gamma-ray hurts (GRBs) - detected by BATSE - into short and long 
subclasses is incomplete. These articles show that the earlier long subclass alone 
should be further separated into a new "intermediate" subclass (2 s < Tg < 10 s) 
and into a "truncated long" subclass (T go > 10 s). Therefore, in what follows, the 
long subclass will contain only the GRBs with T 90 > 10 s, and the intermediate 
subclass will be considered as a new subclass (2 s < T 90 < 10 s). The "short" 
subclass is defined by T 90 < 2 s (for definition of Tg see [3]). 

Fully independently, Balazs et al. [4,5] suggest that GRBs are distributed 
anisotropically on the sky. In addition, they show that the short subclass shows 
an anisotropy, but the intermediate + long subclasses do not show. The different 
behavior is confirmed on 99.7% confidence level [5]. It is difficult to explain such 



behavior of subclasses by the instrumental effects alone. Hence, some intrinsic 
anisotropy should exist. A recent study [6], which is based on the spherical har- 
monic analysis, shows that just the GRBs of "intermediate subclass" have an 
intrinsic anisotropy on confidence level 97%. 

Here we shortly describe and summarize the new results of two further tests. 
The details of these tests will be published elsewhere [7,8]. 

GRBs will be taken between trigger numbers 0105 and 6963 from Current BATSE 
Catalog [3] having defined T 90 (i.e. all GRBs detected up to August 1996 having 
measured T 90 ). From them we exclude, similarly to [9,4,5], the faintest GRBs 
having a peak flux (on 256 ms trigger) smaller than 0.65 photon/(cm 2 s). The 1284 
GRBs obtained in this way define the "all" class. From them there were 339 GRBs 
with T 90 < 2 s (the "short" subclass), 181 GRBs with 2 s < T 90 < 10 s (the 
"intermediate" subclass) and 764 GRBs with T 9 o > 10 s (the "long" subclass). We 
will study the all class and the three subclasses separately. 



We separate the sky in declination into md ec > 1 stripes having the same "effec- 
tive" area (4ir / rridec steradian). This means that the boundaries of stripes are the 
declinations, which ensure that the probability to have a GRB in a stripe is for any 
stripe the same (= l/m dec ). Because the sky-exposure function of BATSE is not 
depending on right ascension, this may easily be done by the convenient choices of 
declinations. We also separate the sky in right ascension a into m ra > 1 stripes. 
Hence, we separated the sky into M = m dec x m ra areas ("cells") having the same 
"effective" size Att/M steradian. 

If there are N GRBs on the sky, then n = N/M is the mean of GRBs at a cell. Let 
Hi] i — 1, 2, M be the observed number of GRBs at the ith cell (Y^fti n i — AO- 
Then 



defines the observed variance. For the given cell structure with M cells, the 
measured variance varu should be identical to the theoretically expected value 
nil - l/M). 

Q cell structures may be probed for the same sample of GRBs. We will choose 
mdec = 2, 3, . . . , 8 and m ra = 2, 3, . . . , 16. I.e. it will be Q = 105. In the coordinate 
system with axes x — l/M versus y = ^jvaru/n = (var /mean) 1 / 2 the Q values of 

(var /mean) 1 / 2 define Q points, where j = 1, 2, Q). The theoretical expectation 
is verified by least squares estimation ( [10], Chapt. 5.3.1.). Our estimator is the 
dispersion 



COUNTS-IN-CELLS TEST 
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var M = (M-l)- 1 Y / (n t -n) 2 



(1) 



i=l 




(2) 



We throw 1000-times randomly iV points on the sphere, and repeat the above 
calculation leading to ctq for every simulated sample. Then we compare the size 
of the oq obtained from this simulation with oq obtained from the actual GRB 
positions. Let u be the number of simulations, when the obtained oq is bigger than 
the actual value of cfq. Then (100 — u/10) is the confidence level in percentage. 



TWO-POINT ANGULAR CORRELATION FUNCTION 

TEST 



Be given No GRBs on the sky. There are Nd(Nd — l)/2 angular distances 
among them. These measured distances are binned into bins with binwidth 4 
degree. Then, using a Monte Carlo simulation, there are distributed randomly N R 
points on the sky, where N R 3> N D . (In our case: N R = 1000 N D .) Then the 
Nr(N r — l)/2 angular distances are binned in the same manner. In addition, the 
NrNd random-data pairs are also binned in the same manner. 

In [11] the following formula is proposed in order to obtain the w{6) correlation 
function: 

< DD> < RD> 
w(6) = — — 2 — — — + 1, 3 

y J <RR> <RR> ' K ' 



where < ... > means a normalized mean (see [11]). 

The la uncertainty in w{9) for the given bin is given by 5w(9) = n iD ' D [11]. This 
formula allows to test the zero expectation value of w(9). For any bin be calculated 
the dimensionless value x from the relation \w(0) — x5w(6)\ = 0. Then xa is the 
probability that w = holds for the given bin. Because w = must occur for any 
bin, simply the biggest value of xa defines the "Poissonian" confidence level [11]. 

We will verify the confidence levels by Monte Carlo simulations, too. We will take 
- instead of the the actual BATSE positions - randomly generated N D positions, and 
then we will repeat the above procedure. This simulation will be done 500 times. 
If there are u such simulations, when the absolute value of w S i m is bigger than the 
absolute value of w act , then (500 — ou)/5 is the confidence level in percentage. 

As the final confidence level the smaller value will be taken. 

We exclude the non-uniform sky exposure function of BATSE instrument as 
follows. Assume that the Monte Carlo simulation defines a point at a given position. 
We generate for any such point an additional random number between and 1, 
too. If this number is bigger than the actual value of sky-exposure function at that 
point, then this point is omitted from the Monte Carlo sample. 

Note that this test is usually more effective than the counts-in-cells one, because 
usually the "distance-based" tests are more sensitive than "cell-based" tests ( [10], 
Chapt.2.6). 



THE RESULTS 



The counts-in-cells tests give u = 287 (uj = 80, uo = 36, oo = 440) for all (short, 
intermediate, long) GRBs. Hence, the rejection of null-hypothesis of randomness 
is confirmed for the intermediate subclass on the 96.4% confidence level. For the 
short and long subclasses, respectively, and also for all GRBs the null-hypothesis 
cannot be rejected on the > 95% confidence level. 

The calculated four correlation functions give the following results: 

1. For the short subclass there is an essential departure from zero for 9 = (14 ±2) 
degrees. We have a 99.2% confidence level for the non-randomness. 

2. In the case of intermediate subclass the "suspicious" angles are the values 
9 = (6 ± 2), 9 = (50 ± 2), and 6 = (90 ± 2) degrees. For 9 = (6 ± 2) degrees we 
have a 99.8% confidence level for non-randomness. 

3. In the case of long subclass for the angle 9 = (94 ±2) degrees we have a 99.0% 
confidence level for non-randomness. 

4. In the case of all GRBs the following angles are "suspicious": 9 = (50 ± 2), 
9 = (94 ± 2), 9 = (130 ± 2) and 9 = (150 ± 2) degrees. For 9 = (94 ± 2) we have a 
99.8% confidence level for non-randomness. 

FIRST CONCLUSION: The intrinsic non-randomness is confirmed on the 
confidence level > 95% for the intermediate subclass by both methods; the an- 
gular correlation function (counts-in-cells) method gives a 99.8% (96.4%) confidence 
level. This supports the results of [6]. 

SECOND CONCLUSION: For the remaining two sub-classes and for all 
GRBs the null-hypotheses of intrinsic randomnesses are rejected only by the angu- 
lar correlation function method on the confidence levels between 99% and 99.9%. 
The counts-in-cells test, similarly to [6], did not reject the null-hypothesis of ran- 
domness for these cases. 

The peak flux = 2 photons/(cm 2 s) (on 0.256s trigger) is practically identical to 
the medium. Therefore, we consider the GRBs having smaller (bigger) peak flux 
than 2 photons/ (cm 2 s) as the "dim" ("bright") subclass ("half) of the intermediate 
subclass. There are 92 GRBs (89 GRBs) at the "dim half ("bright half). 

The 105 "var/mean" tests for these two parts give the results that the "dim half" 
has an intrinsic non-randomness on the 99.3% confidence level; the "bright half 
can well be random. The sky distribution of 92 dim GRBs is shown on Figure 1. 

THIRD CONCLUSION: The intrinsic non-randomness of the "dim half 
of the intermediate subclass of GRBs is confirmed on the 99.3% confidence 
level by the counts-in-cells method. 
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FIGURE 1. Sky distribution of 82 "dim" GRBs of the intermediate subclass in equatorial 
coordinates. 



